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Abstract 



■ We consider modified theories of gravity with a direct coupling between matter and 

£Nj , geometry, denoted by an arbitrary function in terms of the Ricci scalar. Due to such a 

coupling, the matter stress tensor is no longer conserved and there is an energy transfer 
between the two components. By solving the conservation equation, we argue that the 
matter system should gain energy in this interaction, as demanded by the second law of 
thermodynamics. In a cosmological setting, we show that although this kind of interaction 
may account for cosmic acceleration, this latter together with direction of the energy 
. transfer constrain the coupling function. 
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c<3 ■ 1 Introduction 

Cosmological observations on expansion history of the universe indicate that the universe is in 
a phase of accelerated expansion. This phenomenon may be interpreted as evidence either for 
existence of some exotic matter components or for modification of the gravitational theory. In 
the first route of interpretation one can take a mysterious cosmic fluid with sufficiently large 
and negative pressure, dubbed dark energy. In the second route, however, one attributes the 
accelerating expansion to a modification of general relativity. A particular class of models 
that has recently drawn a significant amount of attention is the so-called f(R) gravity models 
(for a review see, e.g., [1] and references therein). These models propose a modification of 
Einstein-Hilbert action so that the scalar curvature is replaced by some nonlinear function 
f{R). Over the past few years, these theories have provided a number of interesting results 
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on cosmological scales. In particular, there exist viable f(R) models that can satisfy both 
background cosmological constraints and stability conditions [2] . Among these cosmologically 
viable models there are some ones which also satisfy solar system constraints under a chameleon 
mechanism [3]. 

In this context, it is recently shown that introducing an explicit coupling between the Ricci 
scalar and matter Lagrangian may explain the flatness of the rotation curves of galaxies [4]. 
Thus, one generalizes the f(R) gravity models as 

S = J d A x^-g{ l -h{R) + [1 + \f 2 (R)]L m } (1) 

where fi(R) and / 2 (-R) are arbitrary functions of the Ricci scalar R and L m is the Lagrangian 
density corresponding to matter systems. The parameter A characterizes the strength of the 
non- minimal coupling of /^(-R) with matter Lagrangian. When A = 0, there is no such an 
anomalous gravitational coupling of matter systems. In this case, the choice fi(R) = 2kR 
with k = {\QitG)~ l gives the standard Einstein-Hilbert action while a nonlinear fi(R) function 
corresponds to the usual f(R) modified Gravity. 

Varying the action with respect to the metric g^ u yields the field equations, given by, 
(f[(R) + 2\&(R)L m )R lu , - \g,MR) = (V„V„ - g^a)(f[(R) + 2\f 2 (R)L m ) 

+ [l + A/ 2 (i2)]7£ (2) 

where the prime represents the derivative with respect to the scalar curvature. The matter 
energy-momentum tensor is defined as 

rpm _ -2 5{y/=gL m ) 

Due to the explicit coupling of matter systems with Ricci scalar, the stress tensor T™ u is not 
divergence free. This can be seen by applying the Bianchi identities V M G^ = to (2), which 
leads to 

l + A/ 2 (i?)' 

The coupling between matter systems and the higher derivative curvature terms describes 
transferring energy and momentum between matter and geometry beyond the usual one already 
existed in curved spaces. Moreover, it can also lead to deviations from geodesic motion in the 
theory described by (1). Recently, there have been some attempts to use such an anomalous 
coupling to address dark matter problem [4] and the cosmological constant problem [5]. Here, 
particular emphasis is made upon those features of the model (1) that could yield accelerated 
expansion of the universe. In fact, there have been already some works on this issue [6] [7]. 
However, in these works the important role of the non-conservation equation (4) is missed and 
the effect of the i?-matter coupling in producing non-linear extra terms in the gravitational 
field equations is studied which leads to accelerated expansion under certain conditions. In 
the present work, particular emphasis is placed on the role of (4) by solving this equation 
in section 2 and taking into account the evolution of matter energy density. By choosing a 



V %™ = u i - T'P^R (4) 
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power law form for the nonlinear function f2(R), we will show that the non-conservation of 
matter energy density means energy transfer between matter and geometry with a constant 
rate. Thus the energy transfer is constrained by the second law of thermodynamics so that 
the latter only allows injection of energy into the matter system. In section 3, we study this 
issue in a cosmological setting. We will show that this one-way energy transfer constrains the 
allowed range of the parameters of the model to produce cosmic speed-up. In section 4, the 
conclusions are presented. 



2 Conservation Law 

As it is clear from (4), details of the energy exchange between matter and geometry depends 
on the explicit form of the matter Lagrangian density L rn . Here we consider a perfect fluid 
energy-momentum tensor as a matter system 

Tfj,v = (Pm + Pm) u ^U u + p m gfiu (5) 

where p m and p m are energy density and pressure, respectively. The four-velocity of the fluid 
is denoted by u^. 

There are different choices for the perfect fluid Lagrangian density which all of them leads to 
the same energy-momentum tensor and field equations in the context of general relativity [8] 
[9]. The two Lagrangian densities that have been widely used in the literature are L m = p m 
and L m = —p m [4] [7] [10] [11]. For a perfect fluid that does not couple explicitly to the 
curvature (i.e., for A = 0), the two Lagrangian densities L m = p m and L m = —p m are perfectly 
equivalent, as discussed in [10] [11]. However, in the model presented here the expression of 
L m enters explicitly the field equations and all results strongly depend on the choice of L m . 
In fact, it is shown that there is a strong debate about equivalency of different expressions of 
the Lagrangian density of a coupled perfect fluid (A ^ 0) [13]. Here, contrary to [10], we will 
take L m = p m as the Lagrangian density of the matter fluid. 

We project (4) onto the direction of the four- velocity which satisfies the conditions = — 1 
and u u V^u u = 0. We also assume that p m = up m with u being a constant equation of state 
parameter. Then, contracting (4) with gives the conservation equation 

Xf (R) 

u^^pm + (w + l)p m VX = - l\'A L m + Pm )u v V v R (6) 

1 + Ahyn) 

We use Friedmann-Robertson- Walker metric given by the line element 

ds 2 = -dt 2 + a 2 (t)( - ^ -. + dn 2 ) (7) 
VI — kr 2 

where a(t) is the scale factor. Homogeneity and isotropy of the universe implies that = 
(1,0,0,0) and r} = T^q = = H where H — ^ is the Hubble parameter and an overdot 
indicates differentiation with respect to the cosmic time t. The expression (6) is then reduced 
to 

p m + 3H(u + l) Prn = — Xf l[ R ) m (Lm + Pm)R (8) 
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It is evident that the fluid energy is not conserved due to the explicit fluid-curvature coupling. 
The right hand side of (8) acts as a source term describing the energy transfer per unit time 
and per unit volume. This is, however, a general statement and there are some situations that 
in spite of such a coupling the right hand side of (8) vanishes. These situations are as follows : 

1. When one simply chooses L m = —p m . We will see below that in this case even though the 
energy is conserved and there is no matter creation (or annihilation), the fluid particles do not 
follow the geodesies of the background metric. 

2. When R = or R = constant. For instance, during inflation in which the scale factor 
exponentially increases the Ricci scalar remains constant. Thus the R- matter coupling does 
not lead to matter creation (or annihilation) in the inflationary phase for any fz(R) function 
and any choice of L m . 

3. For the choice L m = p m and taking p m + p m = we have again vanishing of the right 
hand side of (8) for any fz(R) function. This equation of state belongs to a perfect fluid 
which describes a cosmological constant (with equation of state parameter oo = —1). This is 
important since there is a tendency in the literature to model a cosmological vacuum decay 
scenario by considering an interaction between vacuum and cold dark matter [14]. Thus R- 
matter coupling models can not provide such vacuum decay scenarios. 

We now project (4) onto the direction normal to the four-velocity by the use of the projection 
operator h^ u = u^u u + g^ u . This results in 

h^VT^ = (u + l)p m u v V v u a + V> m + u'VV^p™ 



A/ 2 (i2) 



(L m -p m )h av W u R (9) 



This is equivalent to 



with 



u „w = -— + r^uV = r (io) 



This is an additional force exerted on a fluid element implying a non-geodesic motion. Notice 
that since h au u a = 0, we have f a u a = and the additional force is orthogonal to the four- 
velocity. This is consistent with the usual interpretation of the four-force, according to which 
only the component of the force orthogonal to the particles four-velocity can influence their 
trajectory. 

The additional force due to i?-matter coupling should be attributed to the first term. The 
second term proportional to the pressure gradient does not exhibit a new effect and is the 
usual term that appears in equations of motion of a relativistic fluid. In our choice, L m = p m , 
the first term on the right hand side of (11) vanishes implying that fluid elements follow 
geodesies of the background metric and there is no additional force. In this case, matter is 
still non-conserved and the equation (8) takes the form 

p m + 3H(p m + p m ) = — + l) Pm R (12) 
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To make a closer look at this equation, we assume a power-law expansion for the scale factor 
a(t) = aot m and we adopt /^(-R) = ®R n with a, n, ao and m being constant parameters. 
Putting the latter into (12), gives 

XnaRJ 1 ^ 1 . . • 
p m + 3#(p m + p m ) = - — ^ Rn [u + l)p m R (13) 

In the following, we consider two different cases: 



2.1 the case XaR n » 1 
In this case (13) reduces to 

R 

p m + 3H(p m + p m ) = -n(uj + l )pm^ ( 14 ) 

we have 

H = mf -1 

R = Q(H + 2H 2 ) = Qm(2m - l)r 2 (15) 

R _ 2 H 
R m 

By substituting these results into (14), we obtain the relation 

p m + 3^Hp m = (16) 

where 7 = (1 — ^)(w+ 1). This is a simple differential equation with an immediate solution 
of the form 

Pm = POO" 37 (17) 

where po is an integration constant. Alternatively, this solution can be written as 

Pm = p a- 3(1 ^ )+£ (18) 

with e — ^(uj + 1). This states that when e > matter is created and energy is constantly 
injecting into the matter so that the latter will dilute more slowly compared to its standard 
evolution p m oc a~ 3 ^ +1 \ Similarly, when e < the reverse is true, namely that matter is 
annihilated and direction of the energy transfer is outside of the matter system so that the 
rate of the dilution is faster than the standard one. It is important to note that in an expanding 
universe (m > 0) and for a matter system satisfying weak energy condition (uu + 1 > 0), the 
sign of e or direction of the energy transfer is solely given by n. It is shown [15] that all models 
which investigate possible interactions in the dark sector, the second law of thermodynamics 
requires that the overall energy transfer should go from dark energy to dark matter. This 
means that so long as the curvature is amenable to a fluid description with a well defined 
temperature, it should suffer energy reduction during expansion of the universe if the second 
law of thermodynamics is to be fulfilled. One immediate implication of this argument is that 
the second law of thermodynamics is consistent with n>0forw + l>0 and n < for 
u + 1 < 0. 
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2.2 the case XaR n « 1 
In this case (13) takes the form 



R 

p m + 3H(p m + p m ) = -nx(u + l)p m — (19) 

where x = \aR n . Combining this with (15) gives 

2n 

p rn + 3 #(i _ — x ){ Pm + p m ) = (20) 
3m 

Since x << 1, when |^ remains of order of unity, we have (1 — ^x) ~ 1. Thus 

p m + 3H{p m + Pm ) « (21) 

which gives evolution of matter energy density as the standard one 

p m ~ Poa-^ +l) (22) 

with p being an integration constant. In this case matter is conserved and there is no creation 
or annihilation. 

Before closing this section, we would like to comment on the two above cases. In general, 
the i?-matter coupling implies violation of equivalence principle so that one should keep A 
sufficiently small to ensure that the model satisfies local gravity constraints. One should ac- 
tually tune a\ to reduce the effects of such violation below current experimental accuracy. 
In our case, the choice L m = p m make the extra force attributed to the _R-matter coupling 
vanish and there will be no deviation from geodesies motion. In other terms, test particles 
with different compositions follow geodesies of the background metric. In this case, there 
is no constraint on (aX)^ 1 coming from local experiments and the two cases XaR n » 1 
and \aR n « 1 can be interpreted as two regimes in which the curvature R n is, respec- 
tively, large and small with respect to (aX)' 1 . For (aX)^ 1 ~ 1, since R decreases in an 
expanding universe, the two regimes XaR n » 1 and XaR n << 1 correspond to early 
and late times for n > 0. For n < 0, the reverse is true. We will return to this issue 
later. 



3 Accelerating expansion 

We can recast the equations (2) in a way that the higher order corrections are written as an 
energy-momentum tensor of geometrical origin describing an effective source term on the right 
hand side of the standard Einstein field equations, namely, 

G^v = T™ + T^ v (23) 

where ^ 

T ™» = f[ + 2XtiL m T ^ (24) 
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+(V M V, - 9llv U){f[ + 2\f' 2 L m ) + \f 2 T, u } (25) 
For the metric (7) and in a spatially flat case k = 0, the field equations become 

3H 2 = Pm + p c (26) 

H + H 2 = ~[p m + p c + 3{p m + p c )] (27) 
where energy density and pressure corresponding to curvature are 

! -{- l -{h-f[R) + \ujf 2 Rp m 



f{ + 2\cof^p m " 2 
- ZH{f{R + 2\uftR Pm + 2\ujf 2 p m ) + \f 2 p m } (28) 

^ = f>,o! f/ + 2Aw /2V m ) + (/i + 2Xuf^p m )(R + 3tfi?) 

/i + 2Xuf 2 p m 

- \uf 2 [Rp m - 2(p m + 3Hp m )] + A\uf 2 'Rp m + l -{h - Rf[) + Aw/ 2 p m } (29) 

in which we have set L m = p m = up m . In order to get more realization about the effects of 
nonlinear terms arising from i?-matter coupling, we write the two expressions (28) and (29) in 
two different cases. In one case, they are written for f 2 (R) = 0, 

p c = ^{ 1 -(f 1 -f[R)+3Hf^R} (30) 

Pc = UR 2 f7 + fl{R + 3HR) + hf x - Rf[)} (31) 
J i 1 

In the other case, we consider them for fi(R) = 2kR and f 2 (R) ^ 0, 

A 



2k + 2\uf 2 p r 



-{uf' 2 Rp m - 6uH(f 2 Rp m + f 2 p m ) + f 2 p m } (32) 



Pc = o ZTT 7, i 2R2 f2Pm + V 2 Prn{R + 3HR) - f 2 [Rp m - 2(p m + 3Hp m )] 
2k + 2\u} 2 p m 

+AftRp m + f 2Pm } (33) 

Like usual f(R) gravity models, the former set is written in terms of R and its derivatives 
while the latter set has also terms containing energy density and pressure of matter due to the 
i?-matter coupling. In both sets, p c and p c can be interpreted as energy density and pressure of 
an effective fluid which provides new possibilities in a cosmological setting. A significant part 
of the motivation for f(R) gravity is that it can lead to accelerated expansion at late times 
without the need for dark energy and also at early times without recourse to an inflaton field. 
In fact, under certain conditions which should be met by the function fi(R) in (30) and (31), 
the curvature fluid can take a sufficiently negative pressure and produce a cosmic speed-up. 
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Similarly, the curvature fluid presented in non-minimal coupling models may have a relevant 
role in addressing some problems such as dark matter and dark energy. There is also a hope to 
achieve an explanation for the coincidence problem, due to the appearance of energy density 
and pressure of matter in (32) and (33) [12]. 

For solving (26) and (27), we should first fix the function fi(R). In order to make our analysis 
less complicated and since we are only interested in the effects of i?-matter coupling, we will 
take fi(R) as linear and set fi(R) = 2kR. Moreover, we should have the scaling of p m (t) which 
is given by (17) and (22) for XaR n » 1 and XaR n << 1, respectively. 
For XaR n » 1, the matter system is not conserved and p m is given by 

p m = Poao 37 r 37m (34) 

Taking into account the expressions /^(-R) = oiR n and a(t) = a t m , we obtain 

6m\ap aQ 3l A 

Pc ~ 2K[Qm(2m - l)] 1 -"* 2 ***!™ + 2\anujp a^t 2 

6Xaup a 3y B 

Pr = q 36 

2K[6m{2m-l)] 1 - n t 2n t^ m + 2Xanu p ao 3j t 2 

where 

A = {nw[m(37 + 2) + (2n - 3)] + (2m - 1)} (37) 

4 

B = {-n(n-l)(n-2)+m(2m-l) + 2n(n-l)+mn[-(2m + 2n-3)+7(4n-3) + 3m7(7-l)]} 

(38) 

In the gravitational equations (26) and (27), the left hand side decays as t 2 while time 
variations of the right hand side are given by (34), (35) and (36). Thus, in a curvature 
dominated regime in which p c » p m , one can write 

2 

2n + 3m7 = 2 ^ m = — (1 - n) (39) 

37 

or, equivalently, 

2(nu + 1) 

3(w + l) v ; 

Note that for a dust matter system {oj = 0), this solution gives m — | implying that there is 
no accelerating expansion. 

For XaR n « 1, the matter system is conserved and p m follows the standard evolution p m = 
p cr 3(w+1 ) = p a 3 ( tJ+1 ^-3 m ( w + 1 ). i n this case, p c and p c become 

6mXapoa 3 ^ +l "> A 



2/t[6m(2m - i)]i-«t 2n t 3m ( w+1 ) + 2Xanup ao 3 ^ +1) t 2 
6Xaup a- 3{ul+1) B 

PC = Q/ I -1 \ (42) 

2K[Qm{2m-l)] 1 - n t 2n t^+ l ) +2Xanu poao {l " + ) t 2 
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where 

A = {nw(3mu + 5m - 2) + 2n(n - 1)oj + (2m - 1)} (43) 

4 

5 = {_ n ( n _l)( n _2)-m(2m-l)(n-2)-2n(n-l)(m-l)+mn(w + l)[4(n-l) + (3mw + l)]} 

(44) 

Putting these into (26) and (27) for p c » p m , leads to 

m = 777 7T ( 45 ) 

3(w + l) v ; 

Inserting back (40) and (45) into the equations (26) and (27) gives expressions relating the 
parameters po, ao, po, A and a. 

Accelerating expansion of the universe puts constraints on the parameters oj and n. To see 
this, let us first consider XaR n » 1 which corresponds to early (late) times for n > (n < 0) 
in expansion history of the universe. In this regime, the matter system is not conserved and 
there is an energy transfer between matter and geometry. For m > 1, one can write 

3 1 

As previously stated, the second law of thermodynamics requires that n > (for oj + 1 > 0) 
which translates to— 1 < oo < — 1/3, implying violation of strong energy condition. In Einstein 
gravity, this is the same condition that a perfect fluid (or dark energy) should satisfy in 
order to produce accelerating expansion of the universe. It is also possible that n < (for 
u + 1 < 0). In this case, (46) requires < oj < 1/3 which is not consistent with the second 
law of thermodynamics. 

On the other hand, in the regime XaR n << 1 which corresponds to late (early) times for n > 
(n < 0) the matter is nearly conserved and there is no constraint coming from the second law 
of thermodynamics. For m > 1, the relation (45) gives 

n<-i(3a; + l) (47) 

implying that accelerating expansion is possible for both oj < — 1/3 and oj > —1/3. 
Our power law solutions give a constant equation of state parameter ou c = p c / p c which can be 
written in terms of n and oj. The corresponding parameters space, which is constrained by 
the above conditions for accelerating expansion of the universe, can be subjected to additional 
constraints coming from recent observations on the equation of state of dark energy. To do this, 
we consider XaR n << 1 when oj > —1/3 since it is only in this case that the present model can 
lead to cosmic expansion in the presence of a matter system satisfying energy conditions. One 
can write ou c = uB/mA which for a given oj and using (43), (44) and (45) gives oj c only in terms 
of n. By combining the result with the constraint oj = — 1.02±Q;}g, coming from observations 
on SNe la [16], one can constrain the parameter n. As an illustration, we have plotted oj c in 
terms of n in fig.l for some values of oj. The figure shows that oj c can be in the observational 
bound for n < —5 or n > 6 when < oj < —1/3. In fig.2, oj c is plotted for different values of 
the parameters n and u. In the regions of the parameters space which correspond to oj > , 
oj c can be in the observational bound only when the absolute value of n is large. For n < 0, 
oj c crosses the boundary oj c = — 1 and the curvature fluid can appear as a phantom. 
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4 Conclusions 



In this work we have studied a class of generalized f(R) gravity models in which there is an 
explicit coupling between the Ricci scalar and Lagrangian density of matter systems via the 
arbitrary function /^(-R)- In general, due to this R- matter coupling, the matter stress tensor 
does not remain conserved. Assuming a power-law form for the scale factor and the function 
f2{R), we have solved the (non-) conservation equation in the two cases XaR n » 1 and 
XaR n << 1. In the first case, there is a constant rate of energy transfer from curvature to the 
matter, as required by the second law of thermodynamics. In the second case, however, there 
is nearly no energy transfer between the two components and matter stress tensor is conserved. 
In both cases there is no extra force in the geodesic equation as the choice L m = p m leads to 
vanishing of the first term on the right hand side of the equation (11). 

We then apply the model to a cosmological setting. There are two different cases according 
to evolution of matter energy density. When there is an interaction between matter and 
geometry, the evolution is given by (17) and the exponent n can be both positive and negative 
for u + 1 > and cu + 1 < 0, respectively, as inferred by the second law of thermodynamics. 
However, we have shown that accelerated expansion is possible for — 1 < uj < — | which is 
exactly the same domain for which the cosmic speed-up can be realized for f'2{R) = 0. Thus 
in this case the model does not provide any improvement with respect to Einstein gravity. 
On the other hand, when p m (t) scales as (22) and the matter system is nearly conserved 
there is no constraint on the sign of n coming from the second law of thermodynamics. Thus 
accelerated expansion is possible both for uj < — | and u > — |. The latter case in which there 
is a cosmic acceleration despite the matter part satisfies energy conditions should be regarded 
as the only improvement that the present model provides with respect to Einstein gravity. In 
this case, the equation of state parameter of the curvature fluid is constrained by observations 
so that there is a bound on the parameter n for any given uj. In particular, we have shown 
that the model can not be consistent with observations for u > when the absolute value of 
n takes values of order unity. 
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Figure 1: The plot of u> c in terms of n for some values of u. The red, blue and green lines 
correspond to u> = —0.1,-0.2,-0.3, respectively. The narrow horizontal band indicates the 
observational bound on the equation of state parameter of dark energy —1.21 < u> < —0.89 
coming from observations on SN la. The figure indicates that u c lies in the observational bound 
for n < —6 or n > 5. 
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Figure 2: The plot of u c in terms of different values of n and u. When the absolute value of n 
becomes sufficiently large, u c lies within the observational bound. 
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